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ABSTRACT 

We report the first results from a survey for 1665, 1667, and 1720 MHz OH emission 
over a small region of the Outer Galaxy centered at l ~ 105.0°, b ~ +1.0°. This 
sparse, high-sensitivity survey (AT 4 ps AT m b ~ 3.0 — 3.5 mK rms in 0.55 km s -1 
channels), was carried out as a pilot project with the Green Bank Telescope (FWHM 
ss 7.6 / ) on a 3 x 9 grid at 0.5° spacing. The pointings were chosen to correspond with 
those of the existing 12 CO (1-0) CfA survey of the Galaxy done at a similar resolution 
(8.4'). With 2-hr integrations, 1667 MHz OH emission was detected with the GBT at 
> 21 of the 27 survey positions (> 78%), confirming the ubiquity of molecular gas in 
the ISM as traced by this spectral line. With few exceptions, the main OH lines at 
1665 and 1667 MHz appear in the ratio of 5:9 characteristic of LTE at our sensitivity 
levels. No OH absorption features are recorded in the area of the present survey, in 
agreement with the low levels of continuum background emission in this direction. At 
each pointing the OH emission appears in several components extending over a wide 
range of radial velocity and coinciding with well-known features of Galactic structure 
such as the Local Arm and the Perseus Arm. In contrast, little CO emission is seen in the 
survey area; less than half of the > 50 identified OH spectral features show detectable 
CO counterparts at the CfA sensitivity levels, and these are generally relatively faint. 
There are no CO features without corresponding OH emission in our survey. With very 
few exceptions, peaks in the OH profiles coincide with peaks in the GBT Hi spectra 
(obtained concurrently, FWHM 8 .9'), although the converse is not true. We conclude 
that main-line OH emission is a promising tracer for the “dark molecular gas” in the 
Galaxy discovered earlier in Far-IR and gamma-ray emission, although further work is 
needed to establish the quantitative details of the connection. Further aspects of the 
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OH CO relation are revealed in a scatter plot of CO vs. OH line strengths. This 
plot suggests a rough proportionality between the bright envelope of the CO emission 
(when present) and level of the accompanying 1667 MHz OH emission. Finally, we note 
several cases of anomalous OH emission. One survey position shows several narrow OH 
spectral features which are not well correlated with the Hi spectrum; these features 
have been identified with a nearby known OH-IR star. Also, nine neighboring survey 
positions show enhanced emission at 1720 MHz, consistent with earlier observations and 
with models involving extended regions of elevated particle density. 

Subject headings: ISM: molecules - ISM: structure - local ISM - surveys - Galaxy: disc 
- radio lines: ISM 


Introduction 


The 18-cm radio lines of the OH molecule have been widely observed in the ISM since the first 
det ecti on more than 50 years ago in absorption against the bright continuum emission from Cas A 
( Weinreb et .al. 1 963 ). The first major survey for OH in the Galaxy was also done in absorption 
([Goss 19681 ) against more than two dozen well-known thermal and nonthermal radio continuum 
sources in the northern sky. Within a year or two of the d etectio n in absorpt ion, OH emission 
from several discrete sources in the Galaxy was announced ({Weaver et.al. 19651 ). but it was the 
unexpected appearance of intense and narrow spectral lines of “mysterium” near the OH emission 
lines in several of these sources that captured the most attention at the time. We now recognize 
these features as arising from OH molecules in the ISM which are being excited to rotational levels 
above the ground state by localized non-LTE energetic processes, such as intense IR continuum 
from dust near hot young stars or high transient particle densities in supernova-driven shocks. 
Numerous studies of these bright “OH maser” sources in the Galaxy were undertaken in the years 
following the initial discoveries (for a review, see e.g. Reid & Moran 198 il l. 


In contrast to studies of the relatively-bright discrete sources of OH maser emission, work on the 
extended OH emission in the ISM has been hampered by insufficient sensitivity. As summarized 
Allen et.al. ( 2012 ). early “blind” surveys for_normal OH emissioijj] in the main lines of OH 


m 


at 1665 and 1667 MHz were unsuccessful ( Penzia s 1964; .Knapp & Kerr 1 972 ). and by the early 
1970’s it appeared that, with the equipment available at the time, main-line OH emission from the 
extended ISM could be detected only in localized regions of high density in gas and dust (Heiles 
1968). The first major extensive survey of the northern Galactic plane in all four OH lines was 
done by Turner! (|1979l ) over more than a 6-year period with the 140-ft NRAO telescope beginning 
in 197(|1. Turner’s survey was mainly directed at known sources from many different catalogs of 
stellar optical/IR objects or emission/absorption nebulae; only a small fraction of the pointings 
were not directed towards specific targets. About half of the nearly 1800 pointings show detectable 
OH, often as a mix of emission and absorption. The large velocity extent and roughly-LTE intensity 


1 In this context, “normal” (or “thermal”) refers to emission with LTE line ratios of 5/9 in the main lines at 
1665/1667 MHz. However, the excitation temperature of this emission is gen erally subthermal, i.e. below the kinetic 
temperature of the gas in which the OH molecules are embedded (e.g. Goss 19681 ). 

2 That paper also includes in Table 1 a list of related OH surveys carried out in the period 1966-1976. 
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ratios of many main-line absorption features indicates that these arise in the general Galactic ISM 
subject to Galactic rotation and embedded in the extended Galactic nonthermal radio continuum 
emission. The emission features in Turner’s data appear to be a more eclectic mix of narrower, 
brighter lines characteristic of various maser sources. There are perhaps only one or two dozen 
pointings among the 732 emission sources in Turner’s survey for which the two main lines appear 
to be visible and are plausibly in the expected LTE ratio of 1.8, although the profiles themselves 
are too faint to yield definitive results. In fact, we now know that the sensitivity of Turner’s survey 
was too low to reveal main-line OH emission in LTE from the general ISM: Turner gives the 5a 
sensitivity of his survey as 0.18 K in AT 4 , which becomes 240 mK in A T m ), units appropriate for 
extended sources observed with the 140-ft telescope. Howev er, ty pical peak values for OH main-line 
profiles in the envelopes of molecular clouds are 50-100 mK (jWannier et.al. 1 9931 ) in the same units, 
and more recent observations of nearby high-latitude diffu se cir rus clouds show profiles with peak 
values of only 10-20 mK outside of the CO-bright regions (jBarriault et.al. 20101). 


Our understanding of the distribution of OH in the Galaxy has recentl y bee n greatly improved 
by the results of two recent surveys. First, the SPLASH survey (Dawson et.al. 2014) mapped OH 
in an area of A l x A b = 10° x 4° along the plane of the Inner Galaxy in the southern sky with the 
Parkes Telescope, and found that OH was widely detected in all four of the transitions which were 
observed. Because the brightness of the radio continuum emission concentrated in the Galactic 
plane is often close to the OH excitation temperature, the profiles show a mix of emission and 
absorption. Features of Galactic structure were observed, but it was not possible to detect diffuse 
OH beyond the central CO-bright regions of molecular complexes to a typical sensitivity limit of 
A T m b ~ 30 mK (3cr). The detailed observation in the four transitions allowed the authors to detect 
and identify many OH masers, and to evaluate the excitation conditions in the diffuse OH gas. 


The second survey ( Allen et.al. 2012, 20131 ) examined the diffuse 1667 MHz OH emission in 
a small A l x A b = 1.0° x 4° region of the Outer Galaxy with a greater sensitivity of AT m (, ~ 15 
mK (3<r). The radio continuum background from the Galaxy is significantly fainter here. The 
correspondence between features in the diffuse OH emission and the neutral hydrogen at the same 
position was good, but fewer than 10% of the OH emission features were found in the CfA Survey of 
the 12 CO(1-0) emission in the Galaxy ( Dame et.al.ll200ll ). This may be signific ant because of report s 
over that past decade of evidence for “dark gas” in the Galaxy first reported by Grenier et_ah (2005) • 
In its apparent association with H I, the overall morph ology of the OH emission resembled that of 
the dark gas (for a recent summary, see Langer et.al. 2014j), a potentially important connection 
since the original tracers of that dark gas (Far-IR and gamma-ray emission) did not have discrete 
spectral features and hence could not provide distance information based on Doppler velocities and 
Galactic rotation. 

The present survey has been undertaken as a pilot project with the NRAO Green Bank Tele¬ 
scope (GBT) in part to further explore the apparent connection of the extended OH emission to 
the dark gas in the Galaxy, and in part to begin a study of the 3D spatial structure of this com¬ 
ponent of the molecular ISM. Owing to spectrometer limitations, and to the prese nce o f low-level 
interference in the spectra, the previous blind survey carried out at Onsala by Allen et.al. ( 2012 , 
20131 ) was limited to the 1667 MHz OH line and to distances within ~ 2 kpc of the Sun. The 
present GBT survey, centered on a slightly different direction but still towards the Outer Galaxy, 
records both main lines and takes advantage of the superior sensitivity and spectral coverage of the 
GBT receivers as well as the low interference levels of the National Radio Quiet Zone where the 





























- 4 - 


GBT is located. Finally, the exact correspondence of our OH survey pointing positions to those 
of the CfA CO survey, plus the nearly-identical angular resolutions of the two data sets, allows a 
straightforward observational comparison of these two tracers of molecular gas in the Galactic ISM. 


2. Observing Program 
2.1. Survey construction 


Our pilot survey is located in an area close to the Galactic plane in the second quadrant. It is 
a blind survey, not directed at any specific previously-identified astronomical objects such as Hu 
regions, supernova remnants, etc. However, the location has been chosen for several other reasons. 
First, the line of sight traverses several distinct features of Galactic spiral structure including the 
Perseus and Outer arms. Second, the relation between radial velocity and distance from the Sun 
is very nearly linear, allowing unambiguous determinations of distances and sizes from the mean 
velocities and velocity widths of spectral features. Third, the ambient Galactic radio continuum 
emission near 18 cm wavelength is very faint in this direction (< IK) and gen erall y b elow the 
range o f excitation tem peratures observed for OH in emission (~ 4 — 10K; He ilesl (119691 ): iRieu et.al. 
(1976); Crutcher (1979)), thereby optimizing the detectability of faint Ok| 1. Finally, there are areas 
of the sky here where the 12 CO(1-0) emission is very faint or even absent from the CfA survey, 
enhancing the opportunity to learn more about the “CO-dark” component of molecular gas in the 
Galaxy. 

As described in (JTJ any program aimed at a quantitative measurement of OH spectral features 
will require long integration times in order to achieve sufficient sensitivity for parameter estimation 
on profiles with very modest S/N. As an example, for t heir study of OH main-line emission in the 
translucent high-latitude cloud MBM40 with the GBT, Gotten et ah (2012) spent about 1/2 hour 
of integration at each pointing, but even then they did not achieve enough S/N to use the 1665 
MHz profiles. An improvement in S/N of at least a factor of 2 is required to render the 1665 data 
useful, providing additional independent estimates of profile parameters and helping to identify 
any profiles exhibiting anomalous excitation. We conclude that the integration time requirement 
is at least ~ 2 hours per pointing for the main-line OH spectra with the GBT. Without further 
improvements in receiver sensitivity, even longer integration times would be required to measure 
LTE-level emission from the 1612 and 1720 MHz satellite lines, and we have not considered that 
possibility here. However, the GBT receivers are indeed capable of simultaneous recording of one 
or even both of these satellite lines. We have therefore obtained spectra of the 1720 MHz lin e as 
well, since enhanced emission from this line has been observed widely in the Galactic ISM (Turner 
1 982 ). The results of the Onsala survey also showed that the OH emission has an extent similar to 
the Hi both in space and in velocity; hence Hi profiles at the same pointing positions are also of 
interest. Fortunately, the GBT receivers provide that capability, although the 21-cm HI line is so 
bright that integration times of only a few minutes per pointing were sufficient and were therefore 
carried out separately. 


3 More recent work ( Liszt fe Lucasl 1996l i has shown that the OH detected primarily in absorption profiles has 


significantly lower excitation temperatures, close to the 2.8 K cosmic background. However, for better or worse, this 
component contributes very little to the OH emission. The fraction of OH in this state is presently unknown. 
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In the absence of prior information on the typical linear scale sizes of OH features, the grid 
spacing of 0.5° is a compromise among competing factors including sky coverage, integration time, 
and the limited clock time available on the telescope for such programs. An important result of 
the present survey will be information on the spatial str uctu re of OH emission which we detect and 
a comparison of that with the SPLASH survey (Dawson et.al. 20141 1. which is providing structural 
information on the combination of OH emission and absorption at smaller angular scales (~ 15') 
with somewhat lower sensitivity in the Southern Hemisphere (Inner Galaxy). The number of 
pointings is similarly time-limited, but here there are some guidelines available from the Onsala 
results. A reasonable goal is to obtain ~ 100 OH spectral components fo r a statisti cally-meaningful 
comparison with the corresponding 12 CO profiles. The mini-survey by 1 Allen et.al. (2012) consisted 
of 81 profiles, each with several sub-components. However, the separation of profile sub-components 
at each pointing was severely compromised by low S/N; the few high S/N pointings in that mini¬ 
survey show as many as 3 OH components. With sufficient S/N at each pointing, ~ 100 OH 
components could perhaps be obtained from ~30 pointings. 
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this field is ss 640 MJy/sr at l = 103.74°, b = +2.17°. The source nearest to a survey position is 
~ 410 MJy/sr at l = 104.58°, b = +1.35°. Grid lines are at intervals of A l = 1.00°, A b = 0.50°. 


With the previous considerations in mind, the parameters of our pilot survey are as follows: 
the survey covers a 3 x 9 = 27-point grid centered at l = 105.0°, b = +1.0° with a grid spacing 
of 0.5°. The total integration time is 2 hours per pointing for the OH lines and 5 minutes for Hi. 
Figure QJi shows the survey grid positions superposed on a small portion of the CfA 12 CO(1-0) 
survey. The small circles at each grid point show the approximate FWHM angular resolution of all 










































3 data sets discussed here, OH (7.6'), CO (8.4'), and Hi (8.9'). Figure [T}u shows the 100/U IRAS 
emission in the area. 


2.2. Observations 


The observations were made with the Robert C. Byrd Green Bank Telescope (GBT), using the 
Gregorian receiver system operating in the the frequency band 1.15 - 1.73 GHz (“L-band”). The 
main beam efficiency of this instrument is 0.95 in this frequency band, as determined by NRAO 
staff (Maddalena 2009, Figure 1), and hence the antenna and main-beam brightness temperatures 
T m b = Ta/ 0.95 are virtually identical. Results given in this paper use both scales, but the choice 
will be clarified in the text, caption, or axis legend. The receiver is a single-beam, dual-polarization 
instrument, having an effective system temperature in the range 16 - 20 K, depending upon the 
weather and the level of the Galactic background continuum emission. The signal is fed to the 
control room by means of an IF system for which we chose a bandwidth of 12.5 MHz in order to 
minimize the possibility of harmful radio interference. The signal from the IF system was copied 
and directed to two sections of the GBT Spectrometer, one centered at 1666.4 MHz, the other 
at 1720.5 MHz. Because of the anticipated weakness of the OH line at 1612 MHz, and because 
of the greater likelihood of interference at that frequency, we did not observe that transition. 
Each 12.5 MHz spectral window supported both (linear) polarizations, and the sampling of the 
autocorrelation data at 9 levels and subsequent Fourier transformation provided 8192 channels of 
spectral data sampled at a frequency interval of 1.526 kHz with a frequency resolution of 1.846 
kHz. The IF band centered at 1666.4 MHz is sufficiently wide to record both of the “main” lines of 
OH at 1665.4018 and 1667.3590 MHz; these were separated later during the data reduction phase. 

Because of the expectation that there was likely to be widespread OH emission at a low 
level everywhere in the survey area ( Alle n et.al.1 201.2), we chose to use frequency-switching as the 
basis for removing the instrumental signature from the data rather than position-switching. The 
switching was made by moving the center frequency of the 12.5 MHz band by ±2 MHz from the 
nominal center frequency and subtracting the shifted spectrum. The wide IF bandwidth in use 
meant that with this scheme the signal spectrum appears in both switching cycles, but in different 
channel numbers. The spectra from the two cycles were therefore shifted and inverted (“folded”) 
before subtraction so as to average the signal data from both parts of the switching cycle. While 
optimizing the S/N and avoiding the possible subtraction of signal still present in a reference field, 
this technique samples a different part of the IF bandpass during the two parts of the switching 
cycle and hence can leave residual baseline structure in the output spectrum. We explored the use 
of an additional reference position located at a high Galactic latitude, but this required significant 
time and movement of the telescope and little improvement to the baselines was obtained from 
it. In addition, this runs the risk of inadvertently subtracting signals from very local gas. We did 
balance the receivers before each 10 minute scan, and this appeared to offer some improvement in 
spectral baseline stability. Scans were accumulated such that the total integration time for each 
survey position in the OH lines was ~ 2 hours. A separate observation of Hi at each pointing 
position was made; the instrumental configuration was the same as for the OH observations, except 
that the center frequency was set to be 1420.4 MHz, and only one scan of 5 minutes duration was 
required because of the great strength of the HI signal. The main observing program was completed 
in 14 sessions of ~ 5 hours each between September 2013 and January 2014. The pointing of the 
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telescope was checked at the start of each session using standard procedures for the GBT. 


2.3. Data reduction 

The HI observations required no correction for baselines at the levels of interest here. Also in 
general only one of the two polarizations was reduced since the HI line is known to be unpolarized. 
The OH data required more careful treatment, and each polarization in each 10-minute OH scan was 
reviewed for the presence of interference or other instrumental problems. In general, the quality 
of the data was very high, and little editing was required. The data for each polarization was 
separately averaged and archived. The data for the two polarizations were then averaged, a baseline 
was subtracted, the resulting profile was smoothed by two channels with a Gaussian convolution 
(without decimation), and the final result was archived for analysis. The baseline was determined 
with a polynomial fit over the velocity range of —320 to +200 km s _1 (approximately 1900 channels) 
in segments chosen to avoid signal peaks. The order of the polynomial was progressively increased 
until the improvement in the fit was negligible; a value of 5 was typical for the 1667 spectra. The 
resultant rms noise level in the OH spectra is typically 3.0 - 3.5 mK in T m & units. The final spectral 
resolution of 3.05176 kHz corresponds to 0.64, 0.55, and 0.53 km s _1 at the Hi, 1665/1667 MHz 
OH, and the 1720 MHz OH lines, respectively. Note that the channel spacing in the final smoothed 
but undecimated spectra is half the quoted spectral resolution. 


3. Survey Results 


3.1. General features of the Hi, OH, and 12 CO(1-0) line profiles 


Figure [2] shows a mosaic of resu lts from this pilot survey along with the corresponding 12 CO 
spectra from the CfA survey (jDame et.al. 2001). For each position in the 3x9 = 27-point grid we 
present 3 spectra aligned in radial velocity above each other in a single “stacked-plot” panel, with 
Hi at the top, 1667 MHz OH in the middle, and 12 CO(1-0) at the bottom. As indicated on the 
plots, the y-axis amplitude scales are in units of Ta = 0.95 x T m b for the GBT Hi and OH spectra, 
and T m i, for the CfA CO spectra. This figure shows several major features of large-scale Galactic 
structure at this l, b location including the Local Arm, the Perseus Arm, and (only in H i) the Outer 
Arm. As a specific example, consider the three spectra for l = 104.5°, b = +2.0° on the right side 
of Fig. [2)3. The Hi profile (top panel) extends over almost the entire velocity range observed, from 
— 130 < Vlsr % +15 km s _1 , with five distinct peaks near Vlsr ~ —2 km s~ x (the immediate solar 
vicinity), Vlsr ~ —20 km s -1 (the Local Arm), Vlsr ~ —50 km s -1 and Vlsr ~ —65 km s -1 (the 
Perseus Arm), and Vlsr ~ —100 km s _1 (the Outer Arm). The 1667 MHz OH spectrum shows 
features one can identify with the solar vicinity, the Local Arm, and the Perseus Arm. Noteworthy 
is the fact that in general every OH peak has an associated H 1 peak at essentially the same radial 
velocity, although the converse is not true at present levels of OH sensitivity. Cases where OH 
peaks appear without a corresponding HI peak are rare and can often be identified with OH maser 
sources, as we shall describe in 1. '1.5.2 1 Note also that, at this particular pointing, faint CO emission 
features can be identified near 0, -20, and -50 km s -1 . We will discuss the quantitative comparison 
of OH and CO in © 







3.2. Variations in profile structure on scales of 0.5' 


The profiles displayed in Figure [2] show significant changes over the angular scale of 0.5° used 
in the survey sampling. The changes are less pronounced for the features at Vlsr ~ 0 km s” 1 
because even the cores of molecular clouds would be expected to subtend angles of 0.5° if they 
are at distances of 200 pc. The features which arise from molecular clouds in the Perseus Arm 
show much greater variation. As an example, the strong feature seen at approximately -50 km/s at 
position l = 105.0°, b = +1.0° is weaker by a factor of four in three of the four adjoining pixels. It 
continues to be visible at l = 105.0°, b = +0.5°, but is actually stronger at l = 104.5°, b = +0.5°, 
an d a seco nd feature has appeared at —65 km s -1 in the latter spectrum. Recent determinations 
( Choi et.al. 2014) measure the distance to the inner part of the Perseus Arm in this direction to be 
3.2 kpc. The grid spacing in the present OH data therefore corresponds to a linear scale of ~ 28 
pc, sufficient to resolve broad structure in a giant molecular cloud, but too coarse to resolve core 
structures in such clouds. In order to compare the OH molecular clouds with the more familiar CO 
molecular clouds, it will be necessary to map selected OH features on a finer grid, perhaps even 
down to the Nyquist interval of the GBT beam. 


3.3. Gaussian fits to discrete features in the 1667 MHz OH spectra 

Gaussian fits were found for all distinct OH features, although such fits were often a poor 
representation since many of the profile features are too complex to be effectively fitted by even 
multiple Gaussians. However, the Gaussian fits do provide rough values for profile parameters that 
are useful for discussion and for a general comparison with other ISM components. The 1667 MHz 
Gaussian fits have amplitudes ranging from 7 mK to 60 mK, with the brightest recorded profile 
feature peaking at 100 mK. The average amplitude is near 25 mK with a standard deviation of 15 
mK. Profile widths show a wide range from 1 — 10 km s , with an average of ~ 3.9 km s^ 1 and a 
standard deviation of 2.5 km s -1 . We will discuss the line strengths based on direct integration of 
the profiles in H4.ll 


3.4. 1665/1667 OH emission line strength ratios 


The ratio of the main OH lines at 1665 and 1667 MHz is an important indicator for the 
excitation mechanism of OH molecules in the ISM, and may also be used as a measure of optical 
depth. The satellite lines at 1612 and 1720 MHz are quite sensitive to anomalous excitation by Far- 
IR photons and shocks, but the main lines at 1665 and 1667 MHz usually appear in the LTE ratio 
of 5:9 characteristic of optically-thin collisional excitation. Larger main-line ratios (up to 1:1) have 
been interpreted as indicating higher optical dep th (Heiles 1969), but cases of non-LTE behavior of 
the main lines have also been reported (Grutcher 1979!)• However, accurately determining line ratios 
by comparing the velocity-integrated line strengths can be difficult with weak emission lines, as the 
result is sensitive not only to receiver noise but also to uncertainties in the baseline subtraction 
process. We will discuss this problem in more detail in H4.21 Figure [3] illustrates a quantitative 
but in some ways more robust method for making this comparison. In this Figure we present as 
an example the OH spectra recorded at 3 adjacent pointings in the grid located at l = 105.50° 
(left column), l = 105.00° (middle column), and l = 104.50° (right column), but all at the same 
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latitude of b = +0.50°. This presentation is a variation on the “stacked-plot” presentations of 
Figure [2] only here each column shows the 1665 spectrum along the upper row, the 1667 MHz 
spectrum along the middle row, and the scaled difference spectrum assuming low-opacity emission 
lines in LTE (symbolically, 1665 — 1667/1.8) along the bottom row. Note that the scaled difference 
spectra do not show features that echo the main line profiles (either + or —), thereby strongly 
suggesting that the main lines are indeed generally in LTE even though the baseline levels may 
not be rigorously zero over the relevant velocity ranges. This is true even for profiles which show 
significant substructure on small scales in radial velocity, such as the example shown in the right 
column of Figure [21 Those few pointings where the difference spectra clearly show relatively narrow 
features above the noise will be discussed in sections 13.5.21 and 13.5.41 Only one feature of the more 
than 50 in our survey shows evidence for high optical depth, and a second may be enhanced by a 
nearby source of Far-IR emission; these cases will be discussed in £13.5.31 
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Fig. 3.— Main-line OH intensity ratios for three adjacent pointings at b = +0.50° in the observed grid: 
l = 105.50°, (left column), l = 105.00° (middle column); and 104.50° (right column). The spectra shown are: 
1665 MHz (top row), 1667 MHz (middle row), and the scaled profile differences 1665 — 1667/1.8 (bottom 
row). Note that these scaled differences show no discrete features which mimic the original line profiles. 


3.5. Miscellaneous noteworthy spectra 

3.5.1. A deep exposure at l = 105.5°, b = +3.0° 

In an attempt to extend the OH detections to distances as far as possible from the Sun, we 
have made a longer exposure at the pointing l = 105.50°, b = +3.00° (located in the upper left 
corner of the observing grid in Figured]) with the aim of recording OH emission in the Outer Arm 
of the Galaxy near —100 kms -1 . The Hi-OH-CO stacked plot for this pointing (left column in 
Figure [2ji) is reproduced here in the left column of Figure [4] for comparison. The right column of 
Figured] shows all three long-exposure spectra in a 1665-1667-1720 MHz stacked-plot display. The 
total integration time for the spectra in the right column is 403 min, and the final rms noise after 
baseline subtraction is ~ 1.8 mK. The improvement over our usual 2-hr integrations (left column, 
middle panel) is very nearly a factor 2, as expected, indicating that receiver instability is in general 
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Fig. 4.— A long 403-min exposure spectrum (right column of 3 panels) is shown next to the standard 
1-hour Hi-OH-CO stacked-plot survey result at l = 105.50°, b = +3.00° (left column). The long exposure 
data in the right column is shown in a stacked plot of 1665 OH (upper panel), 1667 OH (middle panel; 
compare to the standard 1-hr survey result directly to the left), and 1720 OH (lower panel). Note that this 
latter spectrum clearly shows a feature at 1720 MHz in the Perseus Arm near -65 km s -1 . 


not the limiting factor in our observations. Unfortunately we have not succeeded in detecting the 
Outer Arm of the Galaxy in OH emission; whether we have just been unlucky with this particular 
pointing or whether the Outer Arm is less rich in OH emission remains an open question. However, 
this long exposure reveals interesting features at 1720 MHz (right column, lower panel) near zero 
velocity (local) and near -65 km s -1 (Perseus Arm) which will be discussed in H3.5.41 


3.5.2. An OH-IR star near l = 105.0°, b = +2.5° 


Figure [5] reproduces in the left column the stacked Hi-OH-CO plot of Figure [2k from our 
survey pointing at l = 105.0°, b = +2.5°. The OH spectrum in the center panel of this column is of 
a nature not seen at any other survey pointing, showing bright “spikey” features that often do not 
correspond with Hi peaks. Referring to the 1665-1667-1720 MHz stacked plot in the right column 
of Figure [5l these features are especially bright at the following approximate velocities: -2, -33 and 
-40 km s^ 1 at 1665 MHz, -2, -11, -33, and -41 km s^ 1 at 1667 MHz, and -2 km s -1 at 1720 MHz. In 
contrast to the usual situation, the Hi profile at this position (FigureEl left column, top panel) does 
not show a peak corresponding with the -41 km s -1 OH feature, suggesting an association at least 
of this OH feature with a compact “non-ISM” source. Further confirmation of anomalous features 
in the OH results at this pointing are shown in the three panels in the center column of Figure El 
where the scaled difference spectrum for the main OH lines is shown as the lower panel. Here it is 
now clear that not only does the -41 km s -1 feature show an anomalous OH main-line ratio, the -11 
km s -1 feature is also anomalous. It is natural to think that these two spikes might be associated, 
and indeed they are. An examination of archival material reveals the presence of a known OH-IR 
star at l = 104. 9°, b = +2. 4°; this source is listed in the “Hamburg Database of Circumstellar OH 
Masers” ( [Engels fe Heidmann 2011, www.hs.uni-hamburg.de/maserdb) and is close enough to be 


detected on the flank of the GBT point spread function (0.13° FWHM). The spectrum recorded 
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Fig. 5.— Survey spectra at l = 105.0°, b = +2.5°. In the left column note that narrow OH features are 
observed without corresponding Hi features, something that is not seen at any other survey position. In 
the lower plot of the middle column, the scaled difference of 1665 and 1667 MHz spectra show two residual 
features, suggesting that these features are not in LTE. These anomalous spectra can be explained by the 
presence of the OH-IR star 104.9+2.4 within the telescope field of view. See text H3. 5. 21 for a discussion. 


Fig. 6.— 
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Spectrum of OH 104.9+2.4 at 1667 MHz with the Nangay radio telescope, from Wolak et al 


on this object with the Nangay Radio Telescope at 1667 MHz by Wolak et all (2012j) is reproduced 
in Figure [6l It shows the characteristic “double-horned” profile of an OH-IR star, with narrow 
features at both -11 and -41 km s , and accounts for half of the bright narrow spikes seen in at 
1667 MHz in Figure [5j OH-IR stars are usually stronger at 1612 MHz, but unfortunately we did 
not obtain a spectrum at that frequency during these observations. 


Curiously, the narrow line (FWHM ~ 1.5 km s -1 ) at -2 km s” 1 in the 1667 MHz spectrum of 
Figure [5] (center column) appears to be in LTE as far as the main lines are concerned, although the 
bottom panel in the right column of Figure [5] shows a bright 1720 MHz feature at this velocity. This 
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feature does not seem to be associated with the OH-IR star, but rather is an example of several 
survey pointings with enhanced emission at 1720 MHz which will be presented and discussed further 
in ffhOl 


3.5.3. Main OH line Non-LTE features 


Figure [7] shows scaled difference stacked plots of 1665 and 1667 main-line spectra at l = 105.00°, 
b = +2.00° (left column) and l = 104.50°, b = +1.50° (right column). Each of these two spectra 
shows one feature that appears to deviate slightly from the expected LTE ratio. The first feature 
is near -20 km s -1 in the left column, where the 1665 line is apparently too bright and the scaled 
difference spectrum shown in the lower panel has a positive residual. This could occur if the 
1667 MHz feature is somewhat optically thick; in that case the usual optically-thin 1665/1667 
MHz LTE line ratio of 5/9 trends towards 1/1. Such a p henom e non was first noted in the early 
OH observations of dense dust clouds in the Galaxy by Heile s (119691 ). who used the measured 
deviation of the line ratio from 5:9 to estimate the opacity. However, this procedure assumes that 
the excitation temperature of the two transitions is the same. In a careful st udy o f absorption and 
emission features in and around a number of continuum sources, [Rieu et.al. (119761 ) found that the 
excitation temperatures measured for the 1665 lines appeared to be a degree or two higher than 
those for the 1 667 lin es, a lthough the authors did not consider such differences to be significant at 
the time. Crutcher (1979) specifically addressed this point with new, more sensitive observations of 
several features in two continuum sources and concluded that such differences were indeed real, and 
furthermore could also occur in the opposite sense, with the 1667 excitation temperature exceeding 
that of the 1665 line. The profiles shown in the right panel of Figure [7] may be an example of 
such an anomaly; the spectra at this pointing show a feature near 1.5 km s^ 1 where the 1667 MHz 
emission is slightly stronger than expected for LTE, leaving a negative residual. 


The origin of these small main-line anomalies is unclear. The models computed bv lGuibert, et.al 


(1978) generally show that the excitation temperatures of the main lines are very nearly equal; how¬ 
ever, there is one model involving warm dust (see their Figure 9) when the 1667 excitation temper¬ 
ature can exceed that of the 1665 line. It is possible that both high optical depth (favoring the 1665 
line) and warm dust (favoring the 1667 line) play roles in acc ounting for the main-line anomalies . 
A tantalizing fact in the present case is that the IRS A data (IMiville-Deschenes fc Lagachel 120051 
reveal a 100 fj, source in close proximity to the line of sight for the pointing illustrated in the right 
panel of Figure [71 This pointing appears to have an overly-bright 1667 line; however, we presently 
have no evidence that the IR source is in fact physically close to the OH-emitting gas we have 
observed. 


3.5.f. Pointings with enhanced 1720 MHz emission 

Although a large majority of the 1665 and 1667 spectral features at the positions we observed 
were found to be in the LTE ratio or close to it, there are at least nine pointings with 1720 MHz 
features that appear to be significantly stronger than the LTE ratio of 1 /9 expected when compared 
to the corresponding 1667 MHz emission at those positions. Some examples are shown in Figure [8] 
These positions show emission features up to 7.7 times stronger than expected for LTE conditions, 
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Fig. 7.— Scaled difference spectra (1665 — 1667/1.8) for two positions that show departures in the main 
1665/1667 OH line ratios from the LTE value of 5:9. Left panel: The 1665 MHz feature near -20 kms" 1 
appears to be too bright, and leaves a positive residual in the scaled different spectrum shown in the lower 
row. Right panel: The 1667 feature near 1.5 km s^ 1 appears to be too bright, leaving a negative residual. 
See text 33.5.31 for further discussion. 


with a mean of 4.4 ±1.1. It is noteworthy that these unusually bright 1720 MHz emission features 
in our survey are all found in the velocity range —10 > Vlsr > 0 km s -1 , suggesting that these 
features are very local to the Sun. Furthermore, the features appear to be clustered into a sub- 
region of our survey area; Figure [9] shows the pointings in our survey grid where these relatively 
strong 1720 MHz features are found. Abnormally bright OH 172 0 MHz emission has been known 
since the early days of 18-cm OH observations of the ISM; see Goss et.al.1 (1 9731) for a brief history 
and the results of an important early study. In a seminal paper, Rieu et.al. (119761 ) showed that the 
excitation temperature of local gas along the line of sight to the radio galaxy 3C123 is negative and 
< — 6 K, indicating a weak maser. These authors also provided an early pumping model to explain 
the emission involving far-infrared radia tion and collisions with neutral and charged particles. The 
models were developed further by Elitzun dl976i) who emphasized the important role of particle 
collisions (and hence elevated density) in pumping the 1720 line, and Turner (1982) subsequently 


sugge sted that these features might be good tracers of spiral arms. More recently, Lo cket t et.al. 


(J1999) have shown that this line can be strongly excited by SN-generated shocks propagating into a 
molecular cloud, and can provide tight constraints on physical conditions in the shocked ISM. This 
model has been successf ully applied by several authors to explain enhan ced 1720 MHz OH emission 
in several known SNR ()Hewitt et.al.1120061 . 120081 : iPihlstrom et.al.1120081 ) . These studies emphasize 
the utility of this emission line as an ISM diagnostic, and we will focus more attention on it in our 
future work. 
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Fig. 8.— Three examples of unusually strong 1720 MHz OH emission in our blind survey. The 1720 MHz 
spectrum is shown in the top row, 1667 MHz in the middle row, and the scaled difference (1720) - (1667)/9 
along the bottom. Left column: a residual feature near V ~ 0 km s , Middle column: a feature near V « 0 
km s , but the Perseus Arm feature near -50 km s _1 appears to be in LTE at this sensitivity; Right column: 
our most sensitive 1720 MHz scaled difference spectrum with 2x better S/N showing an enhanced 1720 MHz 
feature near V « —5 km s -1 . See text N3. 5 .41 for further discussion. 


4. Line Strength Determinations 
4.1. Profile analysis 

The line strengths (profile integrals) of all discrete features which could be visually identified 
on the 27 1667 MHz OH spectra in the survey were calculated, along with the corresponding 
line strengths on the 12 CO(1-0) spectra obtained over the same velocity ranges. The numerical 
expression used is S = ET m j, x AH in units of °K * km s^ 1 , where T m ^ is the main-beam brightness 
temperature, AH is the channel spacing of the data, and the sums are done numerically over the 
channels containing measurable signal for each feature in a spectrum. This latter determination 
is made by visually displaying the 1667 MHz OH spectra and interactively choosing the velocity 
ranges over which the spectrum is to be integrated. Thereafter the 12 CO(1-0) spectrum is extracted 
from the CfA data base and integrated over the same velocity ranges. The procedure was repeated 
until all visually-identifiable discrete 1667 MHz OH features at each pointing were processed. 


4.2. Error analysis 

The ler dispersions in the line strength values were estimated as a combination of receiver 
noise and baseline uncertainty. The cumulative error owing to random receiver noise was computed 
in the usual fashion by determining the value of the rnrs noise in °K per independent spectral 
channel, multiplying by the square root of the number of such channels in each profile integral, and 
finally multiplying by the channel spacing in km s _1 . We have done a 2-channel gaussian smoothing 
on all spectra before archiving, including the spectra on which the line integrals are based. This 
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Fig. 9.— Sketch of the survey area with the locations of enhanced 1720 MHz emission near Vlsr ~ 0 shown 
with diamonds. See text section 111.5.41 for additional discussion. 


has reduced the rms noise in the OH spectra down to the level of 3.0 - 3.5 mK which we quoted 
in 1)2.31 of this paper. However, we did not decimate the final OH spectra; hence neighboring 
channels in the archived spectra at intervals of 0.275 km s~ x are not statistically independent; only 
every second channel at intervals of 0.55 km s^ 1 qualifies for that. The formula for calculating 
the receiver noise contribution E no i se to the error in the line strength S reckoned over the velocity 
interval Vi Q — V} u is therefore E®E ge = 0.0035 x ( Vi a — 14i)/0.55 x 0.275 in units of °K * km s -1 . 
By contrast, the CfA 12 CO spectra have not been decimated, hence the appropriate equation is: 
^noise = 0-32 x y/( Vi Q — 14j)/0.65 x 0.65 where 0.65 is the channel interval in km s -1 of the CfA 
CO spectral data. 

At the other extreme, the baseline fitting procedure used on every spectral feature ensures 
that the contribution of the baseline uncertainties to the errors in the line strength is essentially 
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completely correlated over the entire velocity interval covered by that feature. In order to determine 
this error, the following observational approach was adopted. Sections of empty baseline were chosen 
with a number of spectral channels comparable to that covering a typical emission feature. For all 
OH and CO spectra the sections were integrated and divided by the exact number of channels used 
in order to obtain a per-channel value, and the standard deviations of these values were calculated 
for a set of typical baseline determinations. These standard deviation values are 0.0005 °K per 
0.275 km s^ 1 channel for all of the OH spectra, and 0.05 °K per 0.65 km s^ 1 channel for the CfA 
12 CO spectra. The total baseline error is obtained by multiplying this per-channel value by the 
number of channels in that particular profile, and then by the channel spacing in km s -1 . The 
expressions used are: E^ e = 0.0005 x [( Vi Q — 14j)/0.275] x 0.275 for the OH spectra, and E^ e = 
0.05 x [(Via — Vhi)/ 0.65] x 0.65 for the CO spectra. The final rms errors A S rms on the line strengths 
S are calculated as the quadrature sum of the two error contributions, A S rm s = \JE 2 noise + E 2 ase 
for each of the two tracers OH and CO. 


4.3. A comparison of 1667 MHz OH and 12 CO(1-0) line strengths 


A total of 53 features were visually identified as showing significant 1667 MHz OH emission 
over the 27 pointings of the blind survey. The two features associated with the OH-IR star and a 
third feature possibly associated with a Far-IR source were removed from the sample. The scatter 
plot of Figure [TO] shows the line strengths of the 50 remaining OH features plotted along the X- 
axis with their ±1(7 (horizontal) error bars, many of which are too small to be discerned. The 
corresponding 12 CO(1-0) line strengths with their associated (vertical) error bars are plotted along 
the Y-axis. A horizontal dashed line shows the zero point for the CO data. 


The first point to note about the scatter plot of Figure flOl is that every OH point plotted has 
a S/N in excess of 8; the OH data in this plot is essentially error-free. The second point is that the 
scatter plot appears to show two classes of CO emission: the first is roughly correlated with OH, 
and the second clusters around zero. To enhance this separation we have plotted the 30 low S/N 
CO features (those less than 2 a above the zero line) as filled circles. The third point to note about 
Figure [T0l is that an improv ement in the 12 CO(1-0) sensitivity (such as would be available if we had 
used the FCRAO survey of iHever et ah ( 1998 ) instead of the CfA survey of Dame et.al. (2001) in 
the Outer Galaxy) would undoubtedly reduce the size of the error bars on the CO data, but the 
mean values would in general be at or below the upper limits shown with the error bars in Figure 
[TOl Improved CO data is therefore unlikely to move this second class of CO data points with line 
strengths near zero (filled circles) upwards into the first class of points that are more proportional 
to the OH emission (open circles). We conclude that the 12 CO(1-0) emission accompanying the 
1667 MHz OH features we have measured falls into two roughly-equal categories: CO features that 
show an approximate correlation with OH emission, and CO features that are faint or absent even 
when appreciable S/N OH emission is present. This quantifies the impression obtained from the 
appearance of the stacked-plot spectra in Figure [2j 


It is apparent that the 1667 MHz OH data is revealing a substantially larger fraction of 
molecular gas in the ISM than is indicated by the 12 CO(1-0) line, and this gas is in different 
places. Possible reasons for this behavior include a selective CO photodissociation mechanism 
which leaves OH molecules intact in these regions of the ISM (possibly requiring special conditions 
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that may be unrealistic), or insufficient collisional excitation of the CO(1-0) rotational line as could 
occur in regions of low volume density. A more quantitative discussion including estimates of the 
implied H 2 column densities is in preparation for a subsequent publication. 
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Fig. 10.— Scatter plot of OH 1667 MHz emission line strength (X-axis) and the 12 CO(1-0) line strength 
(Y-axis) integrated over the same velocity ranges. The open and filled circles together show all OH features 
with line strengths in excess of 8cr. Points denoting 12 CO(1-0) data for which the CfA CO line strength is 
within 2cr of the CO zero line are shown as filled circles. 


5. Conclusions 


Using the Green Bank Telescope, we have carried out a blind, high-sensitivity survey for OH 
emission over a sparsely-sampled grid of 27 pointings centered at l = 105.0°, b = +1.0°. OH 
emission at 1667 MHz has been detected in 53 separate features distributed over a tot al o f >, 21 


of the 27_pointings. The results confirm and extend the conclusions first reported by Allen et.al. 


(2012, 20131), namely, that main-line OH emission is ubiquitous in the Outer Galaxy, and that 
this emission is often not accompanied by 12 CO(1-0) emission. Additional results from the present 
observing program include: 


• No spectral features are found in absorption, in agreement with the low levels of 18-cm 
continuum background emission in this direction in the Galaxy. 
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• Observations at both 1665 and 1667 MHz show that the main OH lines are generally found 
with the optically-thin, LTE ratio of 5:9 on a detailed, channel-by-channel basis. This would 
be expected if these lines are excited by collisions in a low-density medium. 

• The spectra show several components which are well separated in velocity. These components 
appear to be associated with well-known large-scale features of Galactic Structure including 
the Local Arm and the Perseus Arm. 

• Substantial changes in the structure of the OH profiles can be seen on the 0.5° angular scale 
of this survey. These changes are especially large for features located in the Perseus Arm, 
where the survey grid spacing corresponds to 28 pc. 

• Gaussian fits to OH emission-line features exhibit amplitudes of (T m b) ~ 25 ± 15 mK, and 
line widths of (AH) ~ 3.9 ± 2.5 km s -1 , where the ± values given are indicative of the la 
dispersions in the quoted means. However, many observed features cannot be fitted in a 
satisfactory manner with single gaussians. 

• Consistent with previous work, extended emission is seen from the 1720 MHz satellite line of 
OH. This enhanced emission extends over ~ 1/3 of our survey pointings and appears to form 
a large structure (or structures) on the sky near Vlsr ~ 0 km s -1 . 

• Narrow-band anomalous emission in both main lines of OH is found at a single pointing out of 
the total of 27, and has been identified post facto with a known OH-IR star at l = 104.9°, b = 
+2.4°. 

• Evidence for other non-LTE main-line ratios has been found at only two positions of the 27 
in the survey. The reasons for these enhancements are not clear; in the first case the 1667 
line may be enhanced by a nearby Far-IR source, and in the second case the 1667 line may 
be diminished if the optical depth is sufficiently large. 

• Emission line strengths were determined for 50 distinct 1667 MHz OH spectral features in 
excess of 8 cr, and compared with the corresponding 12 CO(1-0) line strengths calculated from 
the CfA survey data. 20 of these OH features showed detectable CO emission with line 
strengths > 2a in the CfA survey. The line strengths of these relatively-bright CO components 
appear to be roughly correlated with the corresponding OH features. 

• The remaining « 30 bright OH features identified in the present blind survey show no de¬ 
tectable 12 CO(1-0) emission at the 2a level of the CfA survey. While higher-sensitivity CO 
data would undoubtedly turn many of the CO non-detections into measurements, such data 
is not likely to recover the “missing” CO emission. In particular the filled circles in Figure [T0l 
which now cluster near the dashed zero line of CO line strength, are not likely to move up¬ 
wards into the region presently occupied by the open circles, at least not for OH line strengths 
in excess of ~ 0.2 °K * km s . 

The results described here indicate that high-sensitivity observations of the 18-cm emission 
lines from the OH molecule are a promising new tracer for studying large-scale features in the 
Galactic molecular ISM and for revealing the presence of H 2 in amounts and in locations not 
well traced by the 3-mm line of 12 CO(1-0). More sensitive and more extensive blind surveys of 
the Galactic sky in OH emission are likely to expand and consolidate this result. The spatial 
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structure of the OH emission can be measured especially well in the Perseus Arm, and surveys 
carried out with the GBT in the Outer Galaxy at pointing intervals smaller than that used here 
would provide new insight into the scale sizes and physical conditions of this component of the 
molecular ISM. Differences in the excitation mechanisms and optical depths of OH and CO can 
provide new diagnostic tools for studying structural features of molecular clouds. Blind surveys 
which include the satellite lines of OH could identify hitherto unknown OH-IR stars and help us to 
understand the origin of the line enhancements. The appearance of OH emission in locations where 
Hi is abundant but the 12 CO(1-0) line is weak or absent confirms an association of OH with “dark 
molecular gas”. It also suggests a degree of mixing of the neutral atomic and molecular phases in 
the low-density ISM which deserves further scrutiny. Translating the results of observations such as 
those reported here into quantitative measures of H 2 column density requires careful consideration 
of a number of details including OH line excitation temperatures and relative abundances. These 
topics will be addressed in future publications. 
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Fig. 2.— a). Profiles at b = +3.00° for l = 105.50° (left), 105.0° (middle), and 104.50° (right). 
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Fig. 2.— b). As in FigEJr, for b = +2.50°. 
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Fig. 2.— c). As in Fig Hi, for b = +2.00°. 
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Fig. 2.— e). As in Fig Hi, for b 
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Fig. 2.— f). As in Fig Hi, for b = +0.50°. 
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Fig. 2.— g). As in FigEJi, for b = +0.00°. 
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Fig. 2.— h). As in Fig Hi, for = -0.50°. 


LSR Radial Velocity (km/s) 


LSR Radial Velocity (km/s) 

’ 2 CO(1-0) Profile vom Dame et al (2001) 


Galactic Position 
104.50° -0.50° 




R Radial Velocity (km/s) 



GBT - 1420 MHz HI Profile 


GBT - 1420 MHz HI Profile 


GBT - 1667 MHz OH Profile 



GBT - 1667 MHz OH Profile 



,city (km/s) 
GBT - 1667 MHz OH Profile 



31 Velocity (km/s) 

2 CO(l-0) Profile from Dame et al. (200') 


al Velocity (km 
2 CO(1-C) Proale from Dame e: al. (200') 


R Radial Velocity (km/s) 

12 CO(1-0) Profile ‘rom Dame et al. (2001) 


Galactic Position 1 


Galactic Position 


Galactic Position 1 

- 105.50 -1.00 ! 

^ 6 

105.00 -1.00 ! 


- 104.50 -1.00 ! 

— ^'YvtAAv :: ' k/lA ^ 

a 2 

— 1 


-- 


R Radial Velocity (km/s) 


Fig. 2.— i). As in Fig Hi, for b = -1.00°. 
























































